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Peaking consistently in June for nearly eleven years, the annual modulation signal reported by 
DAMA/Nal and DAMA/LIBRA offers strong evidence for the identity of dark matter. DAMA's sig- 
nal strongly suggest that dark matter inelastically scatters into an excited state split by 0(100 keV). 
We propose that DAMA is observing hyperfine transitions of a composite dark matter particle. As 
an example, we consider a meson of a QCD-like sector, built out of constituent fermions whose 
spin-spin interactions break the degeneracy of the ground state. An axially coupled U{1) gauge 
boson that mixes kinetically with hypercharge induces inelastic hyperfine transitions of the meson 
dark matter that can explain the DAMA signal. 



This letter proposes a new class of inelastic dark matter 
(iDM) models that can explain the annual modulation re- 
ported by DAMA/Nal and DAMA/LIBRA [1 . DAMA's 
signal peaks in early June, consistent with dark matter 
scattering, and has remained in phase for nearly eleven 
years. Moreover, the fractional modulation of the signal 
appears anomalously large, and the nuclear recoil spec- 
trum has a peak near Efi ~ 0(30 keV). 

The hypothesis that dark matter scatters inelastically 
off nuclei into a 0(100 keV) excited state elegantly ex- 
plains the salient features of the DAMA signal [2]. IDM 
models predict nuclear recoil spectra with a characteris- 
tic peak and an 0(1) modulation fraction [3]. The large 
dark matter velocity threshold required by inelastic kine- 
matics also implies that heavier nuclei targets like ^^^I 
in DAMA provide enhanced signal sensitivity relative to 
lighter targets such as ^'*Ge in CDMS. 

In composite inelastic dark matter models (CiDM), 
DAMA's observed signal arises from inelastic hyperfine 
transitions of a composite dark matter particle. (For 
other examples of composite dark matter, see [4j.) We 
illustrate this mechanism with a simple model where the 
majority of dark matter is a meson of a strongly coupled 
SU{Nc) gauge sector that confines near A ~ GeV. These 
mesons are comprised of constituent fermions whose hy- 
perfine interactions split the ground state by 0(100 keV). 
When one constituent quark is non-relativistic, a hier- 
archy between the hyperfine scale and the dark matter 
mass follows inevitably from an enhanced spin symme- 
try. The dark matter couples to a new U{1)a' vector 
boson that kinetically mixes with the Standard Model's 
hypercharge [S]. Another version of iDM with kinetic 
mixing is given in [S]. Axial couplings of the U{1)a' 
to the constituent fermions mediate inelastic hyperfine 
transitions that dominate low-energy nuclear scattering. 

The model considered here has two Dirac fermions, h 
and ^'l, transforming in the fundamental representation 



of the SU{Nc) gauge group. The new U{1)a' couples 
axially to '^Ih.l, each of which have equal and opposite 
unit charge. We introduce a charge-2 Higgs 4>, whose 
vacuum expectation value generates a mass for 'ifH,L and 
the A'. The U{1)a' is non-anomalous for this particle 
content. The dark matter candidate is a ^l^h bound 
state, and its stability can be guaranteed by imposing 
a U{1)h~l fiavor symmetry or Z2 under which ^t/f — > 
—"^H- The Lagrangian for this model is 

^ — ^SM + -C* + -^Gaugc + -Cbroak; (1) 
/^Dark Gauge = ~ -Tr G^^ - -i^^, + £^^,"5^,.. (2) 
= ^Llp'fL + yL^L{4>VL+<l)*VB.)'fL 

-Cbrcak = \D^^\^ -\{\4>\'' -vl)\ 

where 'Pl,r are the left and right Dirac projection opera- 
tors, and D^^L,H = {dy, + igDG^±ie'^zA'^)-^ l,h where 
is the SU{Nc) gauge field and ± corresponding to 
and ^'l, respectively. B'"' is the hypercharge field 
strength, F^," is the U{1)a' field strength, and G^i, is 
the field strength for the confining SU{Nc) theory. In 
this model, the dark quarks have a hierarchy of masses 
™L = Ulv^i < A <C mn — yHV4>, and an A' mass in the 
range 100 MeV < uia' < 20 GeV. 

At momenta beneath A, the model consists of mesons 
and baryons built out of 'i'n and 5*^, listed below. Be- 
cause niH ^ rriL, we classify the states by the number of 
^fl- constituents, Nh- Nh can form bound states by 
anti-symmetrizing their color indices. They have binding 
energies 

Eb oc a?(/i*)mH/iV2 (3) 

where at{^*) is the running 't Hooft coupling of the 
strong gauge sector evaluated at the inverse Bohr ra- 
dius of the bound state, /i* ~ atrriH- At distances 



greater than A~^, the color charge of Nh heavy ^ h^s 
can be screened by Nh hght antiquarks to form a 
dark A^j^-meson, or by Nc — Nh 4'l quarks to form a 
dark iV//-baryon. The h and quarks have antisym- 
metrized cofors, so the spins of same-flavor constituents 
must be symmetrized. The resulting range of spins for 
dark mesons and baryons is 



< JnhM < Nh 



1 



':N, - N^ 



H 



< JN^B < \n. 



(4) 



Due to the spin-spin interactions, the lowest-spin config- 
urations with a given Nh are the lowest-energy config- 
urations. In particular, all dark mesons have a spin-0 
ground state. 

Cosmology doniinantly produces Nh = 1, = 0^ 
mesons tt^ = ^l'I's'^h, as will be shown later. In the 
limit rriL < A <C mn, the = 1" vector = ^LJ^l'^H 
is nearly degenerate with Tr^, and is accessible in low- 
energy TTd scattering. Spin-spin interactions generate a 
mass splitting 



A = M, 



Pd 



kA^ 



(5) 



where ~ mn is the ttj mass, and k is an order- 

unity coefhcient that fixes the relation between A and 
A. The TTd and mesons form a multiplet of the 
SU{2)H-spin that rotates '^hs spin. In the limit mn ^ 
A, SU{2)H-spin is an approximate symmetry, guaranteed 
by Lorentz invariance. This is a familiar phenomenon 
in heavy-quark physics and the enhanced symmetry con- 
strains low energy interactions of TTd and pd mesons [7] . 



long as (f> does not significantly mix with the Standard 
Model Higgs boson, elastic dark-matter-nuclues scatter- 
ing induced by (j) exchange is also negligible. 

Using heavy quark effective theory [7], the SU {2)H-spin 
and Lorentz symmetry of the TTd , pd mesons constrain 
the form of their scattering matrix elements. At leading 
order in relative velocity, frci/c — 10"'^, the — >■ Pd 
matrix elements are given by, 

{pd (p', e) I I TTd (p))^ e;, + ^ (p + p')Vj>- , (8) 

where e^, is the polarization of pd, and Cin ^ 1 con- 
trols dipole scattering. We have dropped terms pro- 
portional to q^^ , which have vanishing contraction with 
the conserved Standard Model electromagnetic current. 
The second term in (|8| leads to scattering enhanced by 

/v'^^yA? ~ 10'^ relative to the first term. 

In terms of relativistic effective operators, the low en- 
ergy interactions can be described as: 

£eff = dinM.,^dVd^4 + "-flld^d^pd^F^r (9) 

+ ^9^7rd^c',.7rdF^r + ... 

Elastic transitions mediated by the dd-operator above 
are velocity suppressed relative to the inelastic ones. The 
axial coupling of A' to the fermionic constituents leads to 
a parity constraint on the interactions, forbidding elastic 
operators of the type 



Cforbiddon = 9^ TTd ^ TTd F^f 



(10) 



CIDM SCATTERING 



The low-energy scattering of the dark mesons arises 
after diagonalizing the kinetic mixing terms in ([2]), and 
integrating out the weak interactions. The dark matter 
constituents couple to the Standard Model via 



C 



Int = 



ITLr 



where sg = sin6'w and cg = cos 9^, Jem is the electro- 
magnetic current, J^o is the neutral Z'^ current and Jd 
is the dark sector U{1)a' current. 



■^d 



(7) 



As discussed in [8j , parity-breaking by strong dynamics in 
the dark sector induces elastic charge-radius TTd -I- SM — >• 
TTd -l- SM scattering, which has new consequences for di- 
rect detection. For simplicity, we will discuss the case 
where the strong dynamics preserves parity in the dark 
sector, so that no TTd -I- SM — TTd -t- SM elastic scatter- 
ing is mediated by the t/(I)yi' current in Q to 0(e). So 



that if otherwise allowed would dominate over inelastic 
scattering. 

Using (|8|, the low-energy inelastic cross-section for TTd 
to scatter off a nucleus of mass toat and electric charge 
Z with nuclear recoil energy Eji is 



I'M 



ER\Fii{E]i) 



where 



/cff = m\,/{cinecegdfK 



,(11) 



(12) 



and Fh is the Helm nuclear form factor used in P|. The 
differential scattering rate per unit detector mass is 



dR pqVq f 3 Vrci J., s da 

Ji UInAU^ Jv^i^(En) ^0 dEji 



dE 



(13) 



where po = 0.3 GeV/cm'^ is the local density of dark mat- 
ter, Wmin(£'fl) is the minimum relative velocity required 
to scatter with nuclear recoil energy Ejj, and f{v;Ve) 
is the dark matter velocity distribution function in the 
lab frame accounting for the Earth's variable velocity Vg. 
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Naively cutting off the velocity profile above the galac- 
tic escape velocity Vcsc, we use the Standard Halo Model 
velocity distribution function, 



fiv;ve) oc (^e 



e "0 \e(y^^^-\^-^^\), (14) 



Given the high uncertainty on the halo velocity distri- 
bution and the sensitivity of iDM models to it [TD], we 
marginalize over the velocity parameters vq and Vesc, con- 
strained to be in the range 150 km/s < vq < 350 km/s 
and 480 km/s < Wosc < 650 km/s. 



DIRECT DETECTION 

In DAMA, dark matter dominantly scatters off ^^^I. 
The modulation spectrum and rate of the combined 
1.17 ton-yr exposure [T] constrain the parameters in this 
model. We perform a global fit by marginalizing over 
the 5 unknown parameters: M^^^, A, fcs, vq and i^osc- 
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FIG. 1: Values of fcft, defined in (12 I, and A, defined in 
that fit the DAM A/LIBRA signal and do not super- 
saturate the null searches at 2a for a benchmark point of 
Mt^^ — 70 GeV are shown in dark green (90%C.L.) and light 
green (99%G.L.). 

Fig. [l] shows a benchmark mass of M^r^ — 70 GeV with 
90keV < A < 150 keV and 50 GeV < Us < 400 GeV. 
The 90% and 99% "confidence level" contours are plot- 
ted, corresponding respectively to < Xq + ^-^ ^^'^ 
< Xo + where Xo corresponds to the best fit point. 

We compute the x^ using the 12 half-keVee bins of 
DAMA's modulated signal and reported error bars be- 
tween 2 — 8 keVee, as well as the combined high energy 
bin from 8-12 keVee. This model fits DAMA's reported 
rate and nuclear recoil spectrum remarkably well. We 
also include in the x^ contraints from null searches sum- 
marized in Tablc[j where for each null experiment we take 
as the standard error the 2a Poisson fiuctuation over the 
number of observed events. Moreover, we require the 
predictions for the null experiments not to supersatute 
the observations at the 95%C.L. 



COSMOLOGY 

In the strongly coupled model of this letter, 'i'^ and 
^ L annihilate too efficiently for their thermal abundance 
to account for dark matter. Consequently, an asymmetry 
must be generated for "dark meson number," um oc uh — 
riL. Dark mesogenesis could in principle be tied to the 
Standard Model's baryogenesis. 

The TTd meson is a simple iDM candidate. Here we 
show that a dominant fraction of the dark meson num- 
ber asymmetry resides in ttj mesons rather than exotic 
mesons, baryons, or mesons. 

Exotic mesons and baryons are seeded by Coulomb- like 
bound states of heavy quarks h formed before confine- 
ment, or created by merging of ttj mesons after confine- 
ment. 

Before confinement, the gluon entropy exponentially 
suppresses h bound-state formation down to a tem- 
perature T* ~ E B / ^'O-is / Um) ^ Eb/SO, where s is 
the entropy density of the Universe. For Nc > 4 and 
Mttj ^ O(10'*)A, T* < A so gluon entropy prevents 
bound-state formation with Nh > 2 until confinement. 

At confinement, ^'^ quark- antiquark pairs nucleate to 
screen the color charge of the 'i'n- Confinement pref- 
erentially leads to the formation of Nh = 1 dark me- 
son over higher-iV/f^ dark mesons or baryons. High-A'^^ 
dark meson formation is negligible because the ^>h are 
dilute at the time of confinement. Formation of Nh = 1 
dark baryons is Boltzmann-suppressed for Nc ^ 4 be- 
cause they have {Nc — 2) more ^' l constituents than dark 
mesons, so they are heavier. 

After confinement, heavy-quark binding occurs via 
TTd + Ti'd Ti'd^^^ + G, where G is a glueball and TTd^^^ 
is an Nh — 2 dark meson. For A ^ 0(1 GeV) and 
niH ~ 0(100 GeV), these reactions are endothermic be- 
cause the glueballs have a mass mc = 0(A) > Eb- The 
binding reactions of "^h require large momentum trans- 
fer Pniin \/ 'Tijrd (™G — Eb) ^ A, SO the binding cross- 
section is controlled by perturbative ^>h dynamics. The 
thermally averaged binding cross-section is paramet- 
rically 



(av) 



-{mG-EB)/T 



Nh 



(15) 



which is Boltzmann-suppressed for endothermic binding 
reactions. 

To summarize, h binding is suppressed by entropy 
for T < A, and by the endothermic Boltzmann factor 
for r > A. In fact, Nh 1 dark mesons have spin-0 
ground states and similar scattering properties to TTd, so 
their abundances can be significant. Constraints arise 
only from Nh = Nc dark baryons, with potentially large 
elastic scattering cross-sections. Baryon formation pro- 
ceeds through a sequence of 4* //-binding reactions, so a 
mild suppression of the binding rate at each stage signif- 
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Experiment 


Element 


Reference 


Effective exposure 


Period of run 


Signal Window 


Obs. events 


CDMS '05 


Ge 


•12 


34 kg-d 


2005 Mar 25 - 2005 Aug 8 


10 -100 keV 


1 


CDMS '07 


Ge 


[13 


121.3 kg-d 


2006 Oct 1 - 2007 Jul 1 


10 -100 keV 





CDMS '08 


Ge 


[14 


194.1 kg-d 


2007 Jul 1 - 2008 Sep 1 


10 -100 keV 


2 


XENONIO 


Xe 


H5 


0.3 X 316.4 kg-d 


2006 Oct 6 - 2007 Feb 14 


4.5 -75 keV 


13 


CRESST-II '04 


W 


.16 


0.59 X 0.9 X 20.5 kg-d 


2004 Jan 31 - 2004 Mar 23 


12 -100 keV 


5 


CRESST-II '07 


W 


.17 


0.59 X 0.9 X 48 kg-d 


2007 Mar 27 - 2007 Jul 23 


12 -100 keV 


7 


ZEPLIN-II 


Xe 


.18 


0.5 X 225 kg-d 


2006 May 1 - 2006 Jun 30 


13.9 -55.6 keV 


29 


ZEPLIN-III 


Xe 


L19 


0.5 X 126.7 kg-d 


2008 Feb 27 - 2008 May 20 


10.7 -30.2 keV 


7 


ZEPLIN-III (iDM) 


Xe 


■20 


0.5 X 63.3 kg-d 


2008 Feb 27 - 2008 May 20 


17.5 -78.8 keV 


5 


XENONIOO 


Xe 


21 


161 kg-d 


2009 Oct 20 - 2009 Nov 12 


7.4 -29.1 keV 






TABLE I; Summary of the data from null experiments. 



icantly suppresses baryon formation. This is discussed in 
detail in [llj . 

The nearly degenerate Nh = 1 meson spin states 
and pd are equally populated at high temperatures. If 
the pd decays only through kinetic mixing with hyper- 
charge, the only kinematically allowed decays are to tt^ 
plus photons or neutrinos. These decays have lifetime 
longer than the age of the Universe [22]. Constraints 
from direct detection of pd ^ ""d de-excitation in nu- 
clear scattering imply a fractional number- density bound 

< 10- 



'T-Pd/^-TTd ^ 10~'^|22J. This constraint is endemic to all 
iDM models coupled to the Standard Model only through 
kinetic mixing, and is troublesome if kinetic decoupling 
of dark matter occurs before T ~ 100 keV. However, in 
CiDM models, pd is de-excited through Pd+Pd T^d+^'d 
scattering, with a large cross section set by the size of 



v) 



A- 



For 



1 GcV, this reaction freezes out when 
/n^d ^ 10"^ where r,pi„ < A is 



the dark meson, 
100 GeV and A - 
exp(-A/Tspin) = 
the asymptotic spin temperature. After structure forma- 
tion begins, up-scattering of TTd into pd can occur in dark 
matter halos, but the TTd — TTd scattering is endothermic 
and the rate too small to re-populate the pd state to an 
observable level today. 

Even though the CiDM self-scattering cross-section is 
large enough to depopulate the pd abundance, it is nev- 
ertheless consistent with current bounds on dark-matter 
self-interactions (see Table I of [53]) 



cr fiCm2 / GeV\^ /100GeV\ 
2x10"^ i^^j 



m 



A 



M, 



TTd 



) 



That is safely beneath the strongest present limits of 
a/m ^ 10~^ cm^/g. 



the model proposed in this letter - for example, baryons 
and weakly coupled bound states ("atoms") j24]. Many 
of these generalizations also naturally posses large cross 
sections to de-excite the 0(100 keV) excited state in the 
early Universe [52]. In contrast to weakly interacting iDM 
models, strongly interacting CiDM naturally leads to low 
spin temperatures and avoids de-excitation constraints. 



The recoil spectrum predicted by (11) differs from con 



ventional iDM spectra by a factor oc En, offering a po- 
tential means of discriminating CiDM from iDM through 
direct detection. Moreover, specific CiDM models pre- 
dict the existence of dark matter sub-components that 
arise from the small fraction of TTd mesons that do pro- 
cess into other dark hadrons. Detection of these dark 
matter sub-components provides a striking signature of 
the underlying dynamics of CiDM [TT] . 

The U{1)a' discussed in this letter has purely axial 
vector couplings to the dark quarks. Other charge as- 
signments give rise to qualitatively different direct detec- 
tion signatures, such as admixtures of elastic and inelas- 
tic scatterings. Many of these alternate charge assign- 
ments generate elastic transitions of dark matter with 
distinctive nuclear recoil spectra in addition to the in- 
elastic transition [55] . 

DAMA prefers fos ^ 300 GeV, resuffing in m^' < 30 
GeV, well beneath the energy frontier. Searches at high- 
intensity e+e^ machines for CiDM through A' produc- 
tion of light L matter may confirm or refute this entire 
class of models [5B]. The coming generation of collider 
searches and direct detection experiments will determine 
the true nature of DAMA's signal, possibly leading to 
irrefutable signatures of dark matter, and may unveil a 
whole new sector of physics. 



DISCUSSION 
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